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A  molecular  motor  utilizes  chemical  free  energy  to  generate  a  unidirectional  motion  through  the 
viscous  fluid.  In  many  experimental  settings  and  biological  settings,  a  molecular  motor  is  elastically 
linked  to  a  cargo.  The  stochastic  motion  of  a  molecular  motor-cargo  system  is  governed  by  a  set  of 
Langevin  equations,  each  corresponding  to  an  individual  chemical  occupancy  state.  The  change 
of  chemical  occupancy  state  is  modeled  by  a  continuous  time  discrete  space  Markov  process. 
The  probability  density  of  a  motor-cargo  system  is  governed  by  a  two-dimensional  Fokker-Planck 
equation.  The  operation  of  a  molecular  motor  is  dominated  by  high  viscous  friction  and  large 
thermal  fluctuations  from  surrormding  fluid.  The  instantaneous  velocity  of  a  molecular  motor  is 
highly  stochastic:  the  past  velocity  is  quickly  damped  by  the  viscous  friction  and  the  new  velocity  is 
quickly  excited  by  bombardments  of  surrounding  fluid  molecules.  Thus,  the  theory  for  macroscopic 
motors  should  not  be  applied  directly  to  molecular  motors  without  close  examination.  In  particular, 
a  molecular  motor  behaves  differently  working  against  a  viscous  drag  than  working  against  a 
conservative  force.  The  Stokes  efficiency  was  introduced  to  measure  how  efficiently  a  motor  uses 
chemical  free  energy  to  drive  against  viscous  drag.  For  a  motor  without  cargo,  it  was  proved  that 
the  Stokes  efficiency  is  bounded  by  100%  [H.  Wang  and  G.  Oster,  (2002)].  Here,  we  present  a  proof 
for  the  general  motor-cargo  system. 

Copyright  ©  2008  H.  Wang  and  H.  Zhou.  This  is  an  open  access  article  distributed  under  the 
Creative  Commons  Attribution  License,  which  permits  unrestricted  use,  distribution,  and 
reproduction  in  any  medium,  provided  the  original  work  is  properly  cited. 

1.  Introduction  and  mathematical  formulation 

Molecular  motors  play  a  central  role  in  many  cellular  functions.  For  example,  an  Fq  motor 
utilizes  the  transmembrane  proton  gradient  to  power  the  ATP  synthesis  in  the  Fi  portion  of 
fhe  ATP  synfhase;  a  Kinesin  dimer  hydrolyzes  ATP  to  drive  intracellular  vesicle  transportation; 
and  a  V-ATPase  hydrolyzes  ATP  to  pump  protons  against  a  large  proton  gradient  to  regulate 
intracellular  acidity.  Understanding  the  operating  principles  of  molecular  motors  is  crucial  to 
comprehending  intracellular  protein  transport  and  cell  motility. 
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Abstract  and  Applied  Analysis 


Due  to  the  small  size  of  molecular  motors,  the  inertia  of  fhe  motor  is  negligible.  As  a 
result,  the  motor  operation  is  dominated  by  high  viscous  friction  and  large  thermal  fluctuations 
from  the  fluid  environment  [1].  Because  of  these  properties  of  molecular  motors,  the  results 
for  macroscopic  motors,  in  general,  do  not  necessarily  extend  to  molecular  motors.  In  both 
macroscopic  motors  and  molecular  motors,  a  chemical  reaction  can  be  used  to  generate  a 
conformational  change  and  an  active  force  af  fhe  reaction  site.  The  conformational  change 
along  with  the  active  force  is  then  mechanically  delivered  to  drive  the  motor  motion.  This 
mechanism  of  generafing  directed  motion  is  called  a  power  stroke  motor  [2,  3].  However, 
in  a  molecular  motor,  a  imidirectional  motion  can  be  generated  by  a  completely  different 
mechanism.  For  simplicity,  we  consider  the  case  where  a  motor  is  restricted  to  moving  in  one 
spatial  dimension,  such  as  moving  along  a  polymer  track.  Suppose  we  use  the  free  energy 
released  from  the  chemical  reaction  to  block  thermal  fluctuations  toward  one  direction.  Then, 
the  motor  will  move  toward  the  opposite  direction.  Here  the  unidirectional  motion  is  not 
directly  driven  by  an  active  force  produced  at  the  reaction  site.  Rather  it  is  directly  driven  by 
the  bombardments  of  surrounding  fluid  molecules.  Of  course,  drawing  thermal  energy  from  an 
isofhermal  environment  to  drive  motor  motion  is  not  sustainable  without  a  free  energy  input. 
The  free  energy  for  blocking  the  backward  fluctuations  and  rectifying  the  forward  fluctuations 
comes  from  the  chemical  reaction.  This  mechanism  of  generating  a  unidirectional  motion  is 
called  a  Brownian  ratchet  [4-7] . 

One  of  fhe  main  differences  between  molecular  motors  and  macroscopic  motors  is 
manifested  in  the  issue  of  efficiency.  For  a  macroscopic  motor,  the  efficiency  is  well  defined 
no  matter  it  is  working  against  a  conservative  force  or  working  againsf  a  friction  force.  For 
a  molecular  motor,  we  do  need  to  distinguish  these  two  cases.  When  a  molecular  motor  is 
working  against  a  conservative  force,  fhe  thermodynamic  efficiency  is  well  defined  and  is  fhe 
energy  outpuf  to  the  external  agent  exerting  the  conservative  force  divided  by  the  chemical  free 
energy  consumption  in  the  motor.  When  a  molecular  motor  is  working  against  a  viscous  drag, 
the  situation  is  completely  different.  It  has  no  energy  output  at  all.  One  way  to  define  efficiency 
in  this  case  is  to  proceed  with  the  apparent  energy  output  based  on  the  average  velocity. 
The  efficiency  defined  this  way,  called  Stokes  efficiency,  is  different  from  the  thermodynamic 
efficiency.  First  of  all,  these  two  efficiencies  are  for  two  different  cases:  the  thermodynamic 
efficiency  is  for  a  motor  working  against  a  conservative  force;  and  fhe  Stokes  efficiency  is 
for  a  motor  working  against  a  viscous  drag.  But  this  is  not  the  only  difference.  In  single 
molecule  experiments,  a  motor  can  be  put  to  work  against  a  conservative  force.  In  a  different 
experimental  setup,  the  same  motor  can  be  put  to  work  against  a  viscous  drag.  For  a  molecular 
motor,  both  the  thermodynamic  efficiency  and  the  Stokes  efficiency  can  be  measured  but  in 
two  different  experimental  setups.  Before  we  can  treaf  fhe  Stokes  efficiency  as  a  valid  efficiency 
measurement,  we  need  to  show  that  it  is  bounded  by  100%.  Since  the  apparent  energy  output 
used  in  the  definition  of  fhe  Stokes  efficiency  does  nof  have  a  thermodynamic  meaning,  the 
Stokes  efficiency  being  bounded  by  100%  cannot  be  argued  simply  from  thermodynamical 
point  of  view.  In  this  paper,  we  prove  mathematically  that  the  Stokes  efficiency  is  bounded  by 
100%.  The  proof  is  based  on  the  Fokker-Planck  formulation  of  mofor-cargo  systems.  Below  we 
will  first  describe  this  mathematical  formulation. 

A  molecular  motor,  in  general,  has  many  internal  and  external  degrees  of  freedom.  One 
of  these  degrees  of  freedom  is  associated  with  the  motor's  imidirectional  motion,  the  main 
biological  function  of  the  motor.  For  example,  a  Kinesin  dimer  walks  along  a  microtubule 
filament  toward  the  positive  end  [8,  9] .  There  are  many  levels  of  models  for  molecular  motors. 
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from  simple  kinetic  models  with  a  few  states  to  all  atom  molecular  dynamics.  We  adopt 
a  modeling  approach  of  intermediate  level  in  which  the  unidirectional  motion  is  followed 
explicitly  and  the  effects  of  other  degrees  of  freedom  are  modeled  in  the  mean  field  potential 
affecting  the  unidirectional  motion  [2,  5,  7,  10].  To  introduce  this  modeling  approach  of 
intermediate  level,  we  start  with  the  simple  case  of  a  small  particle  in  a  fluid  environment, 
restricted  to  moving  in  one  spatial  dimension,  and  subject  to  a  static  potential,  (p{x),  where  x 
is  the  coordinate  along  the  spatial  dimension.  This  situation  is  actually  very  close  to  that  of  a 
molecular  motor.  The  main  difference  is  that  a  molecular  motor  is  driven  by  switching  among 
a  set  of  static  potentials,  each  corresponding  to  one  chemical  occupancy  state.  The  particle 
experiences  are  (1)  the  force  derived  from  the  potential,  (2)  the  viscous  drag  force  from  the 
fluid  environment,  and  (3)  the  Brownian  force  also  from  the  fluid  environment.  Both  the  drag 
force  and  the  Brownian  force  arise  from  collisions  of  the  particle  with  the  surrounding  fluid 
molecules.  The  drag  force  is  the  mean  of  this  stochastic  bombarding  force  and  the  Brownian 
force  is  the  fluctuating  part  of  this  bombarding  force.  The  drag  force  on  the  particle  is  always 
opposing  the  motion  and  is  proportional  to  the  velocity:  drag  =  -^u,  where  u  is  the  velocity 
and  ^  is  the  drag  coefficient  of  the  particle.  Here,  we  adopt  the  convention  that  stochastic 
processes  are  denoted  by  boldface  letters.  The  Brownian  force  on  the  particle  has  zero  mean 
and  is  modeled  as  a  Gaussian  white  noise,  which  is  the  formal  derivative  of  the  Weiner  process 
(the  Brownian  motion) .  The  magnitude  of  the  Brownian  force  is  related  to  the  drag  coefficient 
as 

Brownian  force  = 

where  fcg  is  the  Boltzmann  constant,  T  is  the  absolute  temperature  [11],  and  W(t)  is  the  Weiner 
process  [12].  This  is  a  result  of  the  fluctuation-dissipation  theorem  [13-15].  The  stochastic 
motion  of  the  particle  is  governed  by  the  Langevin  equation 

where  m  is  the  mass  of  the  particle.  Equation  (1.2)  has  two  very  different  time  scales:  the  time 
scale  of  the  motor  forgetting  about  its  current  velocity  (the  time  scale  of  inertia),  and  the  time 
scale  of  moving  over  one  period  of  potential  (pix)  (the  time  scale  of  reaction  cycle).  Dividing 
both  sides  of  (1.2)  by  m  and  noticing  that  quantity  to  =  m/^  has  the  dimension  of  time,  we 
rewrite  (1.2)  as 

§  =  =i[u-sW].  (1.3) 

Here  g(t)  is  a  stochastic  process  given  by 

g(,)  =  -D^.V2D^,  (1.4) 

where  D  is  the  diffusion  coefficient  of  the  particle  and  is  related  to  the  drag  coefficient  by  the 
Einstein  relation:  D  =  ksT [1] .  The  solution  of  (1.3)  has  the  form 

u(t)  =g(t)  +exp  ^^^(u(0)  -g(0))  +  ^l^exp  (g(s)  -  g(t))ds. 


(1.5) 
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Note  that  (1.5)  does  not  provide  an  explicit  expression  for  the  solution  of  (1.3)  since  g(f)  is 
unknown.  Neverfheless,  it  is  clear  from  (1.5)  fhaf  fhe  effect  of  the  initial  velocity  u(0)  decays 
exponentially  with  time  scale  to  =  m/^.  For  that  reason,  we  call  fo  the  time  scale  of  inertia.  In 
general,  the  drag  coefficient  is  proportional  to  the  linear  size  of  the  particle  [1],  and  the  mass 
is  proportional  to  the  volume  of  the  particle.  For  a  spherical  particle  of  radius  r,  fhe  mass,  fhe 
drag  coefficient  and  the  time  scale  of  inerfia  are,  respecfively. 


C  =  6  jrijt, 


(1.6) 


where  p  is  the  density  of  the  particle  and  ri  is  the  viscosity  of  fhe  surrounding  fluid.  In  general, 
fhe  fime  scale  of  inerfia  tg  is  proporfional  to  the  square  of  the  linear  size  of  particle  whereas 
fhe  coefficienf  varies  wifh  the  shape  of  particle.  Thus,  for  small  particles,  fhe  fime  scale  of 
inertia  tg  is  very  small.  For  example,  for  a  lafex  bead  of  diameter  1  pm  in  water,  we  have 
tg  =  56  X  10^®  s  =  56  ps.  A  time  period  of  56  ps  is  much  smaller  than  the  typical  time  scales 
of  chemical  reaction  cycles  in  molecular  motors,  which  are  of  the  order  of  milliseconds.  The 
shorf  time  dynamics  of  fhe  particle  is  full  of  fasf  decaying  of  the  instantaneous  velocity  in  time 
scale  of  tg  and  fast  increasing  of  fhe  insfanfaneous  velocify  excifed  by  fhe  stochastic  Brownian 
force.  These  short-fime  dynamical  behaviors  present  significant  difficulties  for  analysis  and 
compulations.  If  is  desirable  to  capture  relatively  long-time  behaviors  of  the  motor  in  the 
time  scales  of  reaction  cycles,  without  resolving  these  short-time  details.  For  that  purpose,  we 
consider  the  mathematical  limit  of  (1.3)  as  time  scale  tg  goes  to  zero.  As  time  scale  fo  goes  to 
zero,  the  limit  of  the  solution  of  (1.3)  satisfies 


dx 

dt 


-D 


keT 


+  Vm 


dW{t) 

dt 


(1.7) 


The  reduction  from  (1.3)  to  (1.7)  in  the  limit  of  small  tg  is  called  the  Einstein-Smoluchowski 
limit  [16]. 

In  a  molecular  motor,  the  potential  is  not  static.  Instead,  the  potential  changes  with  the 
current  chemical  occupancy  state  of  the  motor.  In  this  paper,  we  consider  the  case  where  a 
motor  has  only  one  catalytic  site  (one  reaction  cycle).  The  extension  of  the  analysis  to  the  case 
of  mulfiple  cafalytic  sites  is  straightforward  buf  tedious.  Let  N  be  the  number  of  chemical 
occupancy  stales  of  the  motor  in  consideration.  Let  j  1, 2, . . . ,  N  j  denote  the  set  of  N  occupancy 
states,  and  let  |  (pi  ix),(p2{x), (pN{x)  ]  denote  the  corresponding  set  of  N  periodic  pofentials. 
Here  (pj{x)  is  the  periodic  motor  potential  when  the  motor  is  in  chemical  occupancy  state  j. 
Mathematically,  (pj{x)  is  a  periodic  fimction  of  period  L  where  the  period  L  is  either  one  motor 
step  or  a  multiple  of  motor  steps.  For  a  Kinesin  dimer  walking  on  a  microtubule  filament,  one 
motor  step  is  about  8  nm.  For  the  rotary  motor  of  Fi  ATPase,  one  motor  step  is  one  third  of  one 
revolution  {2jr/3).  The  mechanical  motion  of  the  motor  is  governed  by  Langevin  equation: 


^  ^  dW{t) 

dt  ksT  dt 


(1.8) 


where  S(f)  denotes  the  current  chemical  occupancy  state,  and  (ps(t)ix)  is  the  periodic  motor 
potential  corresponding  to  S(f).  Equation  (1.8)  can  either  be  viewed  as  one  Langevin  equation 
that  varies  with  the  current  occupancy  state  S(f)  or  be  viewed  as  a  collection  of  N  Langevin 
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equations,  each  corresponding  to  an  individual  chemical  occupancy  state.  Equation  (1.8) 
describes  only  the  spatial  motion  of  the  motor.  Equation  (1.8)  is  not  closed  since  it  depends 
on  the  current  occupancy  state  S(t).  To  close  (1.8),  we  need  to  follow  the  time  evolution  of 
the  chemical  occupancy  state.  In  a  transition  from  occupancy  sfate  A  to  occupancy  state  B,  the 
commitment  time  is  generally  much  smaller  than  the  residence  time.  The  commitment  time 
is  the  time  of  the  actual  transition  process.  More  precisely,  the  commitment  time  is  the  time 
period  during  which  the  system  has  left  state  A  but  has  not  yet  arrived  at  state  B.  If  we  define 
a  transifion  sfate  between  state  A  and  state  B,  the  commitment  time  can  also  be  viewed  as 
the  residence  time  of  the  transition  state.  Because  the  commitment  time  of  transifion  is  small  in 
comparison  with  the  residence  time,  we  can  approximate  the  transition  as  instantaneous.  Thus, 
we  model  the  stochastic  evolution  of  the  motor's  chemical  state  as  a  continuous  time  discrete 
space  Markov  process,  also  called  a  jump  process: 

+  o{6t),  j  f  i, 

1  -  (5f)^fci^i(x(f))  +  o{6t),  j  =  i, 

where  ki^j{x{t))  is  the  transition  rate  from  sfate  i  to  state  j.  Notice  that  the  transition  rate 
depends  on  the  motor  position  x(f).  In  other  words,  in  Markov  process  (1.9),  the  evolution 
of  chemical  state  S(f)  is  affected  by  the  motor  position  x(f).  On  the  other  hand,  in  Langevin 
equation  (1.8),  the  evolution  of  motor  position  x(f)  is  affected  by  the  chemical  state  S(f).  Thus, 
Langevin  equation  (1.8)  and  Markov  process  (1.9)  are  coupled,  and  together  they  govern  the 
stochastic  evolution  of  both  the  mechanical  motion  and  the  chemical  reaction  of  the  motor.  In 
most  of  chemical  reacfions,  for  example,  in  the  ATP  hydrolysis  cycle,  the  system  goes  through 
the  set  of  N  occupancy  states  sequentially.  In  this  case,  the  N  occupancy  states  form  a  loop,  and 
the  only  allowed  transitions  are  either  forward  to  the  next  state  or  backward  to  the  previous 
state  along  the  loop.  Here,  by  convention,  we  represent  the  transition  N— >lbyN— >N+1  and 
represent  the  transition  1  — >  NbyN+1  — >  N.  With  this  convention,  we  have  that  fci^y(x(f))  ^  0 
only  when  j  =  i  +  l  or  j  = 

The  derivation  in  [17]  was  for  motor  systems  without  cargos.  In  many  single  molecule 
experiments,  the  motor  is  not  observed/ recorded  directly.  Instead  a  large  latex  bead  is  attached 
to  the  motor,  the  position  of  fhe  bead  is  observed /controlled  with  the  help  of  a  laser  trap  or 
a  force  clamp  [9,  18-20].  Motor-cargo  systems  also  occur  in  biological  settings.  Eor  example, 
a  Kinesin  dimer  walks  on  a  microtubule  filament,  towing  a  vesicle  toward  the  membrane  for 
exporting.  To  accommodate  these  experimental  and  biological  settings,  we  need  to  study  the 
behaviors  of  motor-cargo  systems.  Eigure  1  shows  a  motor-cargo  system  corresponding  to  the 
experimental  setups  in  [9, 19].  A  Kinesin  dimer  walks  on  a  microtubule  toward  the  plus  end 
and  a  latex  bead  is  linked  to  the  Kinesin  dimer  to  visualize  the  motion  of  fhe  Kinesin  dimer  and 
to  exert  a  force  on  the  Kinesin  dimer.  The  force  exerted  on  the  Kinesin  dimer  is  varied  either  by 
changing  the  drag  coefficient  of  the  latex  bead  [19]  or  by  applying  a  conservative  force  on  the 
bead  using  a  laser  trap  [9].  In  the  recently  developed  2D  laser  trap  [20],  a  force  perpendicular 
to  the  direction  of  motor  motion  can  also  be  exerted  on  the  Kinesin  dimer. 

Let  x(f)  be  the  coordinate  of  the  motor  and  y(f)  the  coordinate  of  the  cargo  along  the 
direction  of  motion  as  illustrated  in  Eigure  1.  As  we  discussed  above,  the  motor  is  driven  by 
switching  among  a  set  of  N  potentials,  each  corresponding  to  an  individual  occupancy  state. 


Pr[S(f  +  6t)  =  j  I  S(f)  =  i]  =  ■ 
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Figure  1:  A  motor-cargo  system  in  single  molecule  experiments,  corresponding  to  the  experimental  setup 
in  [9, 19] .  A  Kinesin  dimer  walks  on  a  microtubule  filament  toward  the  plus  end,  and  a  latex  bead  is  linked 
to  the  Kinesin  dimer  to  visualize  the  motion  of  the  Kinesin  dimer  and  to  exert  a  force  on  the  Kinesin  dimer. 
The  latex  bead  is  either  loaded  by  the  force  from  a  laser  trap  or  loaded  by  the  viscous  drag  from  the  fluid 
media. 


The  stochastic  motion  of  the  motor  is  governed  by  an  overdamped  Langevin  equation 


*  - 


+  \/'2dDM 


dW{t) 
dt  ' 


(1.10) 


where  Dm  is  the  diffusion  coefficient  of  the  motor,  and  -E'{\  -y-  R)  is  the  elastic  force  on  the 
motor  exerted  by  the  elastic  link  between  the  motor  and  the  cargo.  Here  E{x-y-R)  is  the  elastic 
potential  and  R  the  rest  length  of  fhe  elastic  link.  The  resf  length  R  can  be  eliminated  using  a 
change  of  variable  y  =  y  +  R.  So  without  loss  of  generality,  we  assume  R  =  0.  The  stochastic 
evolution  of  the  occupancy  state  is  governed  by  jump  process  (1.9).  In  the  motor-cargo  system, 
the  external  loading  force  no  longer  acfs  direcfly  on  the  motor.  The  external  loading  force  acfs 
direcfly  on  the  cargo.  The  stochastic  motion  of  fhe  cargo  is  governed  by 


dt 


Dc 


[f  +  E'{x-y-R)] 
keT 


+ 


dW(f) 

dt  ' 


(1.11) 


where  Dq  is  fhe  diffusion  coefficient  of  the  cargo,  and  /  is  the  external  loading  force  acting  on 
the  cargo.  Note  that  in  the  motor-cargo  system,  the  motor  does  not  directly  feel  fhe  external 
loading  force  on  the  cargo;  the  motor  sees  the  elastic  force  from  the  cargo.  On  the  other  hand, 
the  cargo  does  not  directly  feel  the  internal  motor  force;  fhe  cargo  sees  fhe  elastic  force  from 
the  motor.  All  communications  between  the  motor  and  the  cargo  are  done  via  the  elastic  link. 
In  experiments,  only  average  quantities  can  be  measured  repeatedly  and  reliably.  All  average 
quantities  can  be  calculated  by  following  the  probability  density  of  fhe  motor-cargo  system.  Let 
Pj{x,  y,  t)  be  the  probability  density  that  the  motor  has  at  position  x,  the  cargo  is  at  position  y, 
and  the  chemical  reaction  is  in  occupancy  state  j  at  time  t.  Mathematically,  pj(x,  y,  t)  is  defined 
as 


Pr 

Pj{x,y,t)  =  lim  — 

»0 

6y-^0 


X  <  x{t)  <  X  +  6x 
y  <  y{t)  <y  +  6y' 


6x  6y 


(1.12) 
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The  set  of  probability  densities  {pj{x,  y,  t)  j  is  governed  by  the  two-dimensional  Fokker-Planck 
equation  corresponding  to  Langevin  equations  (1.10)  and  (1.11)  and  jump  process  (1.9)  [16]. 
The  Fokker-Planck  equation  is  based  on  the  conservation  of  probability 

(x) 

In  the  above,  /j  {x,y,t)  is  the  probability  flux  along  the  x-direction  (the  dimension  of  the 
motor  position)  in  occupancy  state  j  and  it  has  the  expression 


{x,  y,  t)  =  -Dm 


1  d‘^jix,y)  ^  dpj\ 
ksT  dx  ^dxj 


(1.14) 


Oj(x,y)  is  the  total  potential  of  the  motor-cargo  system  in  occupancy  state  J,  which  includes 
both  the  internal  motor  potential  caused  by  the  chemical  reaction  and  the  elastic  energy  in  the 
link  connecting  the  motor  and  cargo.  Oj(x,y)  does  not  contain  the  effect  of  the  external  loading 
force,  which  is  accounted  for  separately.  Oy  (x,  y)  has  the  expression 


c5y(x,y)  =  (pj(x)  +  E(x  -  y). 


(1.15) 


(x,  y,  t)  is  the  probability  flux  along  the  y-direction  (the  dimension  of  the  cargo  position) 
in  occupancy  state  j  and  it  has  the  expression 


J^^\x,y,t) 


-Dc 


keT 


■/  + 


d<Dj{x,y) 


Pi 


(1.16) 


Pi+i/2{x,y,t)  is  the  net  probability  flux  densify  (nef  probability  flux  per  imit  area  in  the  x  -  y 
plane)  along  the  direction  of  reaction  from  occupancy  stafe  j  to  occupancy  state  /  -i- 1.  It  has  the 
expression 


Pj+\/2{x,y,t)  =  kj^j+i{x)pj{x,y,t)  -  fcy+i^y(x)py+i(x,y,f).  (1.17) 


Note  that  the  transition  rate  kj^j+\{x)  does  not  depend  explicitly  on  the  cargo  position  (y). 
The  cargo  can  affect  the  chemical  reaction  indirectly  by  changing  the  motor  position.  In 
the  chemical  reaction  governed  by  jump  process  (1.9),  the  transition  rates  ki^j{x)  cannot  be 
arbitrarily  specified.  These  rates  are  constrained  by  detailed  balance-like  condition,  which 
ensures  that  if  the  system  is  brought  to  an  equilibrium,  then  the  equilibrium  solution  is  given 
by  the  Boltzmann  distribution  and  the  probability  flux  vanishes  everywhere  [11].  Specifically, 
suppose  state  A  and  state  B  have  well-defined  free  energy  Ga  and  Gg.  Then  the  transition  rates 
between  these  two  states  must  satisfy 


kA-^B 

ks^A 


=  exp 


/  Ga-Gb 

V  ksT 


(1.18) 


If  follows  that  the  transition  rates  used  in  (1.17)  are  constrained  by 

'(pjix)  -(pj+i{x)\ 


kj^j+i{x)  _  /Oy(x,y)  -Oy+i(x,y) 
kj+i^j{x)  V  ksT 


=  exp 


ksT 


■ 


(1.19) 
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The  boundary  conditions  for  (1.13)  in  the  x-direction  are  pseudoperiodic: 


Pj{x  +  L,y  +  L,t)  =  pj{x,y,t), 

(pj{x  +  L)  =  (pj{x),  (1.20) 

ki^j(x  +  L)  =  ki^j{x). 


In  the  i/-direction,  (1.13)  extends  from  negafive  infinity  to  positive  infinity.  For  fixed  x  and 
large  y,  fhe  elastic  energy  E(x-y)  is  large  and  the  probability  density  pj{x,y,t)  is  exponentially 
small.  Thus,  at  infinity,  we  have 


Pjix,±oo,t)  =  0, 


dpj{x,±oo,t) 

dx 


dpj(x,±oo,t) 


(1.21) 


The  boimdary  conditions  along  the  reaction  direction  are  also  pseudoperiodic: 


pN+j{x,y,t)  =  pj{x,y,t), 

(pN+j{x)  =  <pj{x)  +  AG,  (1.22) 

kN+i^N+jix)  =  ki^jix), 

where  AG  <  0  is  the  free  energy  change  in  one  reacfion  cycle. 


2.  Stokes  efficiency,  previous  result,  and  a  general  proof 

In  [17],  the  Stokes  efficiency  was  proposed  to  measure  how  efficiently  the  motor  utilizes  the 
chemical  free  energy  to  drive  through  the  viscous  fluid.  Lef  (u)  denote  the  average  velocity 
and  let  r  denote  the  chemical  reaction  rate  (average  number  of  reaction  cycles  per  unit  time) 
of  fhe  motor  system  The  Stokes  efficiency  for  a  motor  system  without  cargo  is  defined  as 

^Stokes  ^(_AG)+/(u)' 

where  ^  is  the  drag  coefficient  of  fhe  mofor.  In  a  motor-cargo  system,  the  motor  part  has  drag 
coefficient  =  ksT/ Dm  and  the  cargo  part  has  drag  coefficient  =  ksT / Dq-  Since  both  the 
motor  and  the  cargo  move  through  the  viscous  fluid  with  the  same  average  velocity,  the  Stokes 
efficiency  for  a  motor-cargo  system  is  defined  as 


(^M  +  ^c)(u)^ 
r{-^G)+  f{vi)' 


(2.2) 


The  denominafor  of  fhe  Stokes  efficiency  is  the  total  free  energy  consumed  per  unit  time: 
r(-AG)  is  the  chemical  free  energy  consumed  in  the  chemical  reaction  per  unit  time  and 
/(u)  is  the  free  energy  input  to  the  motor-cargo  system  from  the  external  agent  exerting  the 
conservative  force.  If  fhe  external  conservative  force  is  against  the  motor  motion,  then  /(u)  is 
negative,  which  means  the  free  energy  is  acfually  oufput  to  the  external  agent.  The  numerator 
of  the  Stokes  efficiency,  (^m  +  fc)(u)^,  has  the  dimension  of  energy  per  unit  time.  But  it  does 
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not  have  a  clear  thermodynamic  meaning.  As  a  result,  Jjstokes  <  1  cannot  be  derived  from 
simple  thermodynamic  arguments.  In  [17],  we  proved  ?2stokes  <  1  for  motor  systems  without 
cargos.  Below  we  present  a  proof  for  fhe  general  motor-cargo  system  described  by  Fokker- 
Planck  equation  (1.13)  with  detailed  balance  constraint  (1.19)  and  boimdary  conditions  (1.20), 
(1.21),  and  (1.22).  Specifically,  we  wanf  to  prove 

aM  +  ^c)(u)'<r(-AG)+/(u).  (2.3) 


Since  all  average  quantifies  can  be  determined  from  the  steady-state  solution,  we  consider  the 
steady  state  of  (1.13)  as  follows: 

=  P--*) 

To  facilifafe  the  analysis  below,  we  introduce  a  new  fimction  Hy(x,  y): 

®,(x,y) 

Hj (x, y)  =  +  In  {pj{x,y)).  (2.5) 

Note  that  in  the  steady  state,  everything  is  independent  of  f.  From  the  pseudoperiodic 
boundary  condition  in  the  reaction  direction,  we  obtain  immediately  that 

AG 

Hjv+y(x,y)  =  Hy(x,y)  +  — .  (2.6) 

fx)  in) 

With  the  help  of  function  Hj{x,y),  we  rewrite  the  fluxes  and  as 


/f  ’  =  -Dm 


(y) 


-Dc 


1  dpj 

ksT  dx  ^  dx 

a®, 

ksT  ^  by 


bHj 

ft- I) -Deft 


by  ksT  )  ■ 


The  average  velocify  and  the  chemical  reaction  rate  have  the  expressions: 


(2.7) 

(2.8) 


/-oo  N 

1  1 

(2.9) 

^00  N 

1  1  ^Jj^^dydx, 

Jo  J-a3j=l 

(2.10) 

Pj+i/2dydx. 

J  0  J  -00 

(2.11) 

In  the  above,  we  have  two  expressions  for  the  average  velocity  (u).  The  first  expression  (2.9) 
is  the  average  velocity  of  the  motor,  and  the  second  expression  (2.10)  is  the  average  velocity 
of  the  cargo.  Since  the  motor  and  the  cargo  are  elastically  linked,  they  must  have  the  same 
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average  velocity.  Substituting  flux  (2.7)  into  the  average  velocity  for  the  motor  given  in  (2.9) 
and  applying  the  Cauchy-Schwarz  inequality  yield 


noo  N  r)J-f 

noo  N 

noo  JV  1  r  ^00  JV  /  dH 


/*L  /•oo  N 


(2.12) 


(x) 

Using  the  fact  that  the  total  probability  is  one  and  converting  back  to  flux  ,  we  get 

.L.oo  N  an,- 

Pseudoperiodic  boundary  condition  (1.20)  implies 

^00  ^00 

I  =1  ■ 

x=0 

Integrating  by  parts  the  right-hand  side  of  (2.13)  and  using  (2.14),  we  have 


(2.13) 


(2.14) 


(u)^  <  Dm 


noo  N  QJ) 

-00  j=1 


dx 


-dy  d3 


(2.15) 


Equation  (2.15)  is  an  inequality  for  (u)  based  on  the  average  velocity  of  the  motor.  Now  using 
the  average  velocity  for  the  cargo  given  in  (2.10)  with  flux  (2.8),  applying  the  Cauchy-Schwarz 
inequality,  converting  back  to  flux  and  then  converting  back  to  (u),  we  obtain 


(u)'  = 


•L  POO  N 


0  J  -oo  ;=i 


\  \  dy  ksT 


t2 


dy  dx 


(2.16) 
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Integrating  by  parts  and  using  boundary  condition  (1.21)  lead  to  another  inequality  for  (u): 


^co  N  f 


(2.17) 


Combining  inequalities  (2.15)  and  (2.17)  gives  us 


(u)"  = 


Dr 


Dm  +  Dc 

DcDm 
Dm  +  Dc 


(u)'  + 


D 


M 


Dm  +  Dc 


(u)' 


,  DcDm  f  .  > 
^  DM  +  DckRT^"^^' 


(2.18) 


Up  to  now,  we  have  not  yet  used  the  steady-state  Fokker-Planck  equation  (2.4).  Using  (2.4), 
applying  summation  by  parts,  and  using  relation  (2.6),  we  have 


(u)'< 


-DcD 


cr-'M 


noo  N 

noo 

XP;+i/2(Hy+i  -  Hj)dydx 

-00  j=l 


DcDm  f  (u) 


Dm  +  Dc 
DcDm 


Dm  +  Dc 

J 

DcDm  (-AG) 


(2.19) 


Dm  +  Dc  ksT 

=  Qi  +  Q2  +  Qs- 


nco 

pi/idy  dx  + 

-00 


DcDm  /(u) 


Dm  +  Dc  ksT 


Let  us  examine  the  three  terms  Qi,  Q2,  and  Q3  defined  above  one  by  one.  For  Q3,  we  use  fhe 
relation  {DcDm/ {Dm  +  Dc))-{l/kBT)  =  1/(^m  +  Cc)  to  write  it  as 


Qs 


1 


•/(u). 


Cm  +  Cc 

In  Qi,  converting  the  integral  back  to  the  chemical  reaction  rate  gives  us 

1 


Q2  = 


•r(-AG). 


Cm  +  Cc 

For  Qi,  we  first  use  detailed  balance  (1.19)  to  write  flux  density  Pj+1/2  as 

,  /.  ^i+iW  Fy+i\ 

,  /,  /^;>i(x) -^y(x)  P;+i\\ 

'  <=/-/..«  (1  -  «P  ( — Crf —  yr)> 

Then  we  use  definition  (2.5)  to  write  Hy+i  -  Hj  as 


(2.20) 


(2.21) 


(2.22) 


HjM  -  Hj  = 


<pj+iix)  -  tj>j{x)  ^ 


keT 


Pi 


(2.23) 
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Using  these  two  results,  we  write  Q\  as 

T~)  rL  roc  N 

Qi  =  -  exp  (AHy+1/2))  AHy+i/2di/dx,  (2.24) 

L>M  +  lJcJo 

where  AHy+1/2  =  Hy+i  -  Hj.  Since  the  expression  (1  -  exp(a))a  is  always  nonpositive  for  any 
value  of  a,  it  follows  fhat  Q\  is  always  nonpositive:  Qi  <  0.  Finally,  substituting  Qi,  Q2,  and  Q3 
into  (2.19),  we  arrive  at 


(u)'<^^^^(r(-AG)+/(u))  (2.25) 

which  leads  directly  to  the  desired  result  (2.3).  Thus,  the  Stokes  efficiency  is  indeed  bounded 
by  100%  for  a  motor-cargo  system. 

In  summary,  for  the  general  case  where  a  motor  is  elastically  linked  to  a  cargo,  we 
have  proved  mathematically  that  the  Stokes  efficiency  is  bounded  100%.  Therefore,  fhe  Stokes 
efficiency  is  a  justified  efficiency  measure  for  motor-cargo  systems. 
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such  dynamical  systems  can  exhibit  a  very  important  and  still 
unexplained  phenomenon,  called  as  the  Fermi  acceleration 
phenomenon.  Basically,  the  phenomenon  of  Fermi  accelera¬ 
tion  (FA)  is  a  process  in  which  a  classical  particle  can  acquire 
unbounded  energy  from  collisions  with  a  heavy  moving  wall. 
This  phenomenon  was  originally  proposed  by  Enrico  Fermi 
in  1949  as  a  possible  explanation  of  the  origin  of  the  large 
energies  of  the  cosmic  particles.  ITis  original  model  was 
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